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ON THE DECISION PROBLEMS FOR LINEAR 
RECURRENCE SEQUENCES 

MIN SHA 


Abstract. In this paper, we show that for almost all linear re¬ 
currence sequences of algebraic numbers with integer coefficients 
(or rational coefficients), the Skolem Problem is decidable by only 
testing finitely many terms, whose amount is bounded explicitly 
by using the initial data. The same is also true for both the Posi¬ 
tivity Problem and the Ultimate Positivity Problem of real linear 
recurrence sequences of algebraic numbers with integer coefficients 
(or rational coefficients). 


1. Introduction 

1.1. Background. Linear recurrence sequences appear almost every¬ 
where in mathematics and computer science, and they have been stud¬ 
ied for a very long time; see [11] for a deep and extensive introduction. 
In this paper, we focus on several decision problems relating to such 
sequences over the complex numbers C. 

Recall that a linear recurrence sequence of order m > 1 is a sequence 
{u n }^ = o with elements in C satisfying a recurrence relation 

(1.1) Ui+m ^l^n+m—1 T ‘ ' T CL m U n {jt 0; 1) 2, . . .), 

where aq,..., a m 6 C, a m ^ 0 and Uj ^ 0 for at least one j in the range 
0 < j < m — 1, and the sequence {u n } does not satisfy a relation of 
type (1.1) of smaller length. 

Here, we call Gq,.. . ,a m the coefficients of the sequence {u n }, and 
the initial data of {u n } are aq,..., a m , u 0 , ..., w m -i- Moreover, the se¬ 
quence {u n } is called a linear recurrence sequence of algebraic numbers 
(resp. integers, rational numbers, real numbers) if all its initial data 
are algebraic numbers (resp. integers, rational numbers, real numbers). 

In the sequel, we abbreviate both ‘linear recurrence sequence’ and 
‘linear recurrence sequences’ to ‘LRS’. 
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(1.2) f(X) = X m - ai X m ~ l - a m = l[(X - cq) ei G C[X] 


with distinct aq, a 2 , • ■ •, (which are called the characteristic roots of 
{u n }) and ei > 0 for 1 < i < k. Then, u n can be expressed as 


k 



(1.3) 


i =1 


where /* is some polynomial of degree e* — 1 (i = 1,2,... ,k). We 
call the sequence { u n } simple if k = m (that is e± = ■■■ = e m = 1) 
and non-degenerate if cq/aq is not a root of unity for any i ^ j with 
1 < i,j < k. If {ti n } is non-degenerate, then there are only finitely 
many integers n such that u n = 0. 

For the above LRS {u n }, the three decision problems we study are 
the following: 

1. The Skolem Problem (SP): does u n = 0 for some nl 

2. The Positivity Problem (PP): is u n > 0 for all nl 

3. The Ultimate Positivity Problem (UPP): is u n > 0 for all but 
finitely many nl 

Actually, PP and UPP only make sense for real LRS. Note that SP is 
widely open in general case. 

The above three problems (and related variants) have applications in 
various areas, such as theoretical biology, software verification, proba¬ 
bilistic model checking, quantum computing, discrete linear dynamical 
systems, as well as combinatorics, formal languages, etc; see [21] for 
references. We refer the reader to [20] for a survey on these problems. 

There are only few results towards decidability of SP. For such se¬ 
quences of order 1 and 2, this problem is relatively straightforward. 
Decidability for LRS of algebraic numbers and of orders 3 and 4 is in¬ 
dependently settled positively by Mignotte, Sliorey and Tijdeman [18], 
as well as Vereshchagin [25], while that of order not less than 5 is not 
known. Decidability of SP is also listed as an open problem and dis¬ 
cussed by Tao [24, Section 1.9]; see also [13] for a survey. To taste 
the difficulty of the problem, we want to point out that Blondel and 
Portier [4, Corollary 2.1] showed that it is NP-hard to decide whether 
a given integer LRS has a zero. However, it is known that determining 
whether a given integer LRS has infinitely many zeros is decidable [3]. 

Most recently, Ouaknine and Worrell have made breakthroughs on 
PP and UPP. They have showed that UPP for integer LRS of order 5 or 


3 


less is decidable in polynomial time [21], and PP for simple integer LRS 
of order 9 or less is decidable [22], Moreover, UPP has been shown to 
be decidable for simple LRS of rational numbers and of all orders [23] 
in a non-constructive sense; that is, a given LRS of rational numbers 
can be certified ultimately positive, but no index threshold is provided 
beyond which all terms of the sequence are positive. Furthermore, it 
has been proved that it is NP-hard to decide whether a given integer 
LRS is non-negative; see [1, Theorem 1], 

1.2. Our results. In this paper, we show that the above three prob¬ 
lems can be decidable for almost all LRS of algebraic numbers with 
integer coefficients by testing finitely many initial terms, whose num¬ 
ber can be bounded explicitly by a function which is poly-exponential 
in the initial data. The same is also true for those with rational coef¬ 
ficients; see Theorem 3.5. We also give estimates in some sense on the 
number of such LRS for which SP might be not decidable explicitly; 
see Section 3.5. 


2. Preliminaries 

In this section, we gather some definitions and results which are used 
later on. 

2.1. Basic definitions. Given a polynomial 

/(X) = a 0 X m + chX™- 1 + • • • + a m = a 0 (X - a,) ■ • • (X - a m ) G C[X] 

of degree m > 1 , we assume that the roots (listed with 

multiplicities) are labelled so that |«i| > |cv 2 1 > ••• > |a m |. In case 
|cii| = ••• = |«fc| > |afc+i|, we say that / has exactly k roots with 
maximal modulus. 

In order to interpret the meaning of “almost all” in the results, we 
introduce the following definition. 

Definition 2.1. Given a proposition P related to integer polynomials, 
for integers m > 1 and H > 1, we define the sets 

Sm(H) = {/(X) =X m + ail ™- 1 + • • • + a m G Z[X\ : 

|eij| < H, i — 1,..., m} 

S* m (H ) = {/(X) =a 0 X m + cqX" 1-1 + • • • + a m G Z[X] : 

a 0 7 ^ 0, |eij| < H, 2 = 0,1,..., m}. 


and 
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We say that P is true for almost all monic integer polynomials if for 
any integer m > 1 the following limit holds: 

r \{f e S m (H) : P is true for f}\ 

- = L 


Similarly, P holds for almost all integer polynomials if for any integer 
m > 1 the following limit holds: 


lim 

H^roo 


I {/ g S* m (H) : P is true for f}\ 

s* m m 


1 . 


Since the behaviour of a LRS depends heavily on its characteristic 
polynomial, if a certain property is true for almost all monic integer 
polynomials (resp. integer polynomials), then we say that the cor¬ 
responding property is true for almost all relevant LRS with integer 
coefficients (resp. rational coefficients). Note that every LRS with 
rational coefficients uniquely corresponds to an integer polynomial of 
content 1. 


Definition 2.2. For a polynomial f(X ) e C[X] of degree greater than 
0 , it is called dominant if it has a simple root whose modulus is strictly 
greater than the moduli of its remaining roots. If f(X ) is dominant, 
the root with maximal modulus is called the dominant root of f(X). 
The polynomial f(X) is called non-degenerate if the quotient of its any 
two distinct roots is not a root of unity. 

The following are two results concerning non-degenerate or dominant 
integer polynomials: the first result comes from [7, Theorem 1.1] and 
[8, Theorem 1.1], and the second one is from [7, Theorem 1.4] and [9, 
Theorem 1.1]. 

Lemma 2.3. Almost all monic integer polynomials are non-degenerate 
and dominant. 


Lemma 2.4. Almost all integer polynomials are non-degenerate and 
either have a dominant root, or have exactly two roots with maximal 
modulus (counted with multiplicity). 

Furthermore, several algorithms given in [8] test whether a given 
integer polynomial is dominant or not, and they can be easily amended 
to test whether a given integer polynomial has exactly two roots with 
maximal modulus. 


2.2. Previous results. For SP, Mignotte, Shorey and Tijdeman have 
proved the following general result; see [18, Corollary 1], 
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Lemma 2.5. Let {«„} be a LRS of algebraic numbers defined by (1.1). 
For its characteristic roots in (1.2), without loss of generality suppose 
that let! | > |a 2 | > ••• > \otk\- Define an integer r by 

|Qii| — |a 2 | = • • • = \ot r \ > |a r+ i|. 

If either r = 1, or 2 < r < 3 and at least one of the numbers oti/otj 
with 1 < i < j < r is not a root of unity, then SP is decidable for {u n }. 

Combining Lemma 2.5 with Lemma 2.3, it is easy to see that for 
almost all LRS of algebraic numbers with integer coefficients , SP is 
decidable; from Lemma 2.4, the same is also true for those with rational 
coefficients. Now, the remaining problem is to estimate the computa¬ 
tional complexity. That is, the goal is to obtain a lower bound for the 
index n beyond which every term u n is non-zero. 

The following result is a weak form of [2, Theorem 2]; see also [5, 
Theorem 1] and [12, Theorem 1.3] for the previous results. 

Lemma 2.6. Let {u n } be a real LRS defined by (1.1), and let f(X) 
be its characteristic polynomial. If f does not have a positive root with 
maximal modulus, then there are infinitely many integers n such that 
u n > 0 and also infinitely many integers n such that u n < 0. 

Lemma 2.6 tells us that for PP and UPP, we only need to consider 
those LRS having a positive characteristic root with maximal modulus. 


2.3. Mahler measure and Weil height. Given a polynomial f(x) = 
aox m + • • • + a m = a 0 (x — afi) ■ ■ ■ (x — a m ) G C[x] of degree m > 1, its 
length is defined by 


L(f) — |do| + ''' + 1°: 


its height by 


H(f) = max leiil, 

0 <i<m 


and its Mahler measure by 

m 

M(f) = |a 0 | n max{l, |a«|}. 

i =1 


For each / gC[x] of degree m > 1, these quantities are related by the 
following inequality 


(2.1) 


Hlf)2 < M(f) < H(f)V^+ 1, 



6 


MIN SHA 


for instance, see [26, (3.12)]. We also need the so-called Landau’s 
inequality 

( m 

Y. 

i =0 

which was proved, for instance, in [17]. 

Accordingly, for an algebraic number a, its Mahler measure M(a) 
is defined as the Mahler measure of its minimal polynomial / over the 
integers Z, that is, M(a) = M(f). 

For a number held K, we denote by Mr the set of all valuations v of 
K extending the standard infinite and p-adic valuations of the rational 
numbers Q: \2\ v = 2 if v G Mk is Archimedean, and \p\ v = p _1 if v 
extends the p-adic valuation of Q. In particular, if the valuation v of 
K corresponds to a prime ideal p of K lying above a prime number p, 
we also denote the valuation | |„ by | | p , then for any a G K we have 

H p =p-“ d p(«)/ep j 

where ord p (a) is the exponent of p appearing in the prime decomposi¬ 
tion of the fractional ideal aOx, Ok is the ring of integers of K , and 
e p is the ramification index of p over p. For any v G M K , let K v be the 
completion of K with respect to the valuation v, and let d v = [K v : Q„] 
be the local degree of v. When the valuation v corresponds to a prime 
ideal p lying above a prime number p, we also denote K v by A' p and 
QL by Q p , respectively. 

For the above number held K, the (Weil) absolute logarithmic height 
of any non-zero a G K is dehned by 

(2.3) h(a) = d~ l d v logmax{|a|„, 1}. 

v£M k 



Actually, we have 

(2.4) h(a) = d -1 logM(a), 

see [26, Lemma 3.10] for a proof. 

Given non-zero a G K, in view of (2.3) and h(a) = h(a _1 ), for any 
valuation v G Mk we have 

log|a|„| < dh(a)/d v < dh(a). 


(2.5) 
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In the sequel, we use the following formulas without special reference 
(see, e.g., [26]). For any neZ and /3i, ■ ■ - , (3k, 7 G Q, we have 

h(/?i + • • • + 0k) < h(/3i) + • • • + h(/3fc) + log k, 

h(/3i---/3 fc ) <h(A) + --- + h(/3 fc ), 
h(7”) = MMt), 
h(171) < %), 

h(C) = 0 for any root of unity £ G Q. 

Here, Q stands for the algebraic closure of Q. 

We also need the following result, which is exactly [26, Lemma 3.7]. 

Lemma 2.7. Let f G Z[xi,..., x n ] be a non-zero polynomial in n 
variables. Then, for any algebraic numbers 71 ,..., 7 ™, we have 

n 

Hf( 7l> ■ ■ ■ , 7n)) < lo g^(/) + ll( ^ 7 ^ deg ^ f' 

i —1 

where deg,,, / zs the partial degree of f related to x t . 

2.4. Separating the moduli of roots. The following lemma is a 

classical result due to Cauchy; see [19, Corollary 8.3.2], 


Lemma 2.8. Let f(X) G C[X] be a polynomial of degree m > 1 defined 
by 

f (X) = aoX m + a\X m 1 + • • • + a m , 
where a 0 7 ^ 0 and (ai,..., a m ) 7 ^ (0,..., 0). Then, for an arbitrary 
non-zero root x of f, we have 


mino<i< m | cti \ 

H(f)+ mino<i< m |a* 


< Ixl < 1 + -— r maxllai],..., |a m |}. 
l a o| 


We reproduce [8, Lemma 2.4 and Lemma 2.5] as follows. 


Lemma 2.9. Let f(X) G Z[X\ be a quadratic polynomial. Suppose 
that f has two real roots a and /3 with |a| 7 ^ \(3\. Then, we have 

Wal-tfW^Hif)- 1 . 


Lemma 2.10. Let f(X) G 7L\X\ be a polynomial of degree m > 2, and 
let a and (3 be two roots of f satisfying |a| 7 ^ |/3|. Then, 

(2.6) ||a| - |/3|| > 2 m ( m - 1 )/ 4 ( m+ l)- m3 /4+3rn/4-3 iJ ^)-m3/2+m 2 +m/2-2 

if both a and (3 are non-real. If, furthermore, a is real and (3 is non- 
real, then 

(2.7) 

|| a | _ |^|| > 1 )( m_ 2 )/ 2 ^ m 
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Finally, if both a and f3 are real, then 

(2.8) |H -|/3|| > (2m + 1 )- 3m H(f) 2 ~ 4m . 

Note that for large enough m, (2.8) is better than (2.7), and (2.7) 
is better than (2.6). However, for small integer m, this might be not 
true. For simplicity, we put them together into two uniform forms. 

Lemma 2.11. Let f(X) € 7L\X\ be a polynomial of degree m > 3, and 
let a and (3 be two roots of f satisfying |a| ^ |0|. Then, 

(2.9) 

|| a | _ I^H > 2 ( m + l)~ m3 / 4 + 3m / 4 - 6 //( ^-m 3 +m 2 +m/2-2. 

if furthermore a is real, then 

(2.10) ||a| - |0|| > 2- m ( m - 1 )( m - 2 )/ 2 ( m + i)- m ( m -i)-6 i7 ^)-2m( m -i)-i 

Proof. Note that (2.6) is automatically contained in (2.9). Since for 
m > 4, 

—m 3 /4 + 3m/4 — 6 < — m 2 + 3m/4 — 6 < — m 2 + m — 1/2, 
and for m > 3, 

—m 3 + m 2 + m /2 — 2 < —2m 2 + m /2 — 2 < — 2 m 2 + 2 m — 1 , 

we can see that (2.7) is implied in (2.9) by considering m = 3 individ¬ 
ually. In addition, noticing 

2 _ m(m _l)(m— 2 )( m + 1 ^-m 3 /4+3m/4-6 2" 3m ( m + l)" 3m < (2m + l) _3m 

for m > 4, we know that (2.8) is included in (2.9). 

The inequality (2.10) can be proved similarly by merging (2.7) and 

( 2 . 8 ). □ 

2.5. Bounding the coefficients. For further deductions, we need to 
estimate the coefficients in (1.3) when the sequence {w n } is a simple 
LRS of algebraic numbers. 

Lemma 2.12. Let {u n } be a simple LRS of algebraic numbers of order 
m > 2 defined by (1.1), and let f(X) be its characteristic polynomial. 
Write u n as 

m 

U n ^ ^ bj (\ j , 

3 = 1 

where a\,..., a m are distinct roots of f and all bj are non-zero. Then, 
for any 1 < j < m we have 

m—l g 1 

h (bj) < ^ h(iti) + 2 m ^ h(ct fe ) + m 2 h(<x,) + -m 2 — -m — 1 ; 

7 ^ 2 2 
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furthermore if f e Q[X], let 6 be the smallest positive integer such that 
Sf(X) G Z[X], denote Sf(X) by f*(X) and write 


nx) = j2<-x m -‘ez[x], 


i =0 


then we have 


m— 1 


(2.11) h (bj) < Y h K) + 2 m2 lo S Y I 2 ) + 9 m2 “ 9 m “ L 

4 = 0 \ 4=0 

Finally, if f e Q[X] zs irreducible, we have 

m —1 


(2.12) h(6j) < h(«i) + —m log Y 


* 12 


a,- 


H— m 2, - m — 1. 


4 = 0 


.4=0 


Proof. Here, we follow the arguments in the proof of [10, Theorem 3.1]. 
Notice that 


(2.13) (U 0 , Mi, ... , Mm— i) = (&i, b 2 , . ■ ■, b m ) 


f 1 Oil 

1 «2 

V 1 Qim 


_ 771—1 \ 
«i \ 
_ m— 1 

cr 2 


a. 


771—1 


and ai,..., a m are distinct. To solve the above system of m linear 
equations in m unknowns bi,..., b m , we denote the appearing Vander¬ 
monde matrix by V = (a] _1 ) 1< . m . By [14, Formula (6)], the inverse 
of V is given by V" 1 = (^p) 1 < iJ < m , where 

( f ) 2 ^m—i (cii, . . . , aij, ... , CXm) 

Wij = -- 

J j — 1 m 

n (<*j - «n) n (<Xk-<Xj) 

44=1 k=j +1 


and (jfc («i ,... ,dfj,..., a m ) stands for the h-tli elementary symmetric 
function in the m — 1 variables a i,..., a m without ay; for instance, 
in the case j = m, we have <Ti(ai,..., a m -i) = ay + • • • + a m -i and 
(T m _i(a?i, ..., cx rn —i ) aq • • • 

So, for any j with 1 < j < m we have 

m 

bj ^ ^ 'U'i—l'Wjj . 

2=1 


Since ..., a},..., a m ) is a polynomial with coefficients 1 in 

m — 1 variables ai,..., a m (without aj) of degree m — i, length 





10 


MIN SHA 


and degree 1 in each variable a kl k 7 ^ j, by Lemma 2.7 we find that 

h(a m - i (a 1 , ..., d},..., a m )) < log ( m 1 J + h (a k ). 

\ m 1 ) 


On the other hand, we observe that 
h 11^ ~ “i) ) - “ a J') 


< 


22 (h(«fe) + h(aj) + log 2 ) 


k^j 


Y h(o!fc) + (m — l)h(oij) + (m - 1 ) log 2 . 


Thus, we obtain 


h(t%) < li(afc) + (m — l)h(Q'j) + (m — 1) log 2 + log 
Hence, for 1 < j < m we conclude that 

771 

1 ) + h(uy,)) + logm 

7—1 
777—1 

< 22 h(uj) + 2 m 22 h(cbc) + rn{m — l)li(ay) 

i=0 k^j 

+ 2 m(m — 1) log2 + logm 

m— 1 g y 

< 22 h(wj) + 2m 22 h(afc) + m2 h(%) + a m2 ~ _ 1, 

t=0 k^j 

where we also use the fact that the binomial coefficient ( m_1 ) < 2 m_1 
for any 1 < i < m and 


m — 1 
m — i 


3 1 

2 m(m — 1 ) log 2 + logm < -m 2 — 2 m 


1 . 


This gives the first desired upper bound. 

Now we assume that / is a polynomial over Q. For any 1 < i < m, 
denote by deg a* the degree of a; as an algebraic number, since 


h(Q' i ) 


log M(ati) 
deg a* 
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we deduce that 

1 / m 

h(aji) < log M(ai) < log M(f*) < - log I J^|a *| 2 

\i=0 

where the last inequality follows from Landau’s inequality (2.2). Then, 
collecting relevant results yields the inequality ( 2 . 11 ). 

Finally, the inequality (2.12) comes from the fact that deg a* = m 
for any 1 < i < m when / is an irreducible polynomial over Q. □ 

In fact, slightly better results can be achieved for such sequences of 
order 2 . 



Lemma 2.13. Let {u n } be a simple LRS of algebraic numbers defined 
by (1.1) of order 2. Let f(X) and f*(X) be defined as in Lemma 2.12, 
and write u n as 

u n = feitr” + & 2 «2 5 
where oq,a 2 are roots of f. Then, we have 

2 

h( 6 j) < h(ti 0 ) + h(ui) + | log(£ \a *| 2 ) + | 

i=0 

forj = 1 , 2 . 


Proof. Since for n > 0 


u n = &iaq + 6 2 « 2 > 


we derive 


Mi - u 0 a 2 , , u 0 a i - Mi 

b\ —- ana 02 — - 

ot\ — a 2 — 0^2 

So, we obtain 


h(6i) < 1 i(mi - m 0 «2) + h(«i - a 2 ) 

< 1i(m 0 ) + h(Mi) + h(cti) + 2h(a 2 ) + 2 log 2 

< 1i(m 0 ) + h(Mi) + ^ log(J^ \a *| 2 ) + ^ 

i=0 

Similarly, we have 

2 

h(& 2 ) < h(M 0 ) + h(ui) + ^ log(£ \a*\ 2 ) + ^ 


□ 
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2.6. Linear form in the logarithms of algebraic numbers. One 

key technical tool in this paper is Baker’s inequality on linear form in 
the logarithms of algebraic numbers. Here we restate one of its explicit 
forms due to Matveev [16, Corollary 2.3]. 

First, recall that for a non-zero complex number 2 , the principal 
value of the natural logarithm of 2 is 

log 2 = log | 2 | + \f—l • Arg 2 , 

where Arg 2 is the principal value of the argument of 2 (— 7 r < Arg 2 < 
7r). Note that the definition here coincides with the natural logarithm 
of positive real numbers. We also want to indicate that the identity 
log( 2122 ) = log Z\ + log 22 can fail in our settings. 

Let 

A = bi log + b 2 log a 2 H-b b k log a k , 

where k > 2, b\,...,b k G Z, and or,..., a k are non-zero elements of 
a number field K. Let D = [K : Q] and B = max{|6 i|,..., |&^|}. For 
any 1 < j < k, choose real number Aj such that 

Aj > max{.Dh(aj), | logo'll}. 

Suppose that A ^ 0. Then, we have 

(2.14) log |A| > -2 6k+20 D 2 A 1 ■■■Ak log (eD) log (eB), 

where e is the base of the natural logarithm. 

3. Main results 

3.1. Characteristic polynomials with dominant roots. We first 
want to get explicit versions of some results in Lemma 2.5 in Section 
3.1 and in Section 3.2. 

Theorem 3.1. Let {u n } be a real simple LRS of algebraic numbers 
defined by (1.1) of order m > 3, and let f(X) be its characteristic 
polynomial. Suppose that f(X) is a polynomial over Q and has a dom¬ 
inant root, and denote by d the degree of the number field generated by 
u 0 ,..., u m —i over Q. Let 5 be the smallest positive integer such that 
Sf(X) G 7L\X\. Denote 5f(X) by f*(X), and write 

m 

f(i) = ^<r-’GZ[i]. 

i =0 

Then, for the sequence {u n }, SP, PP and UPP are all decidable by only 
testing the first [W(ii)J terms, where 

Nfiu) = 2 m ( m -P( m - 2 V 2 (rn+l) m{m - 1)+6 H(f*) 2m{m - 1 '> +2 C 1 (u) 
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and 

(m-l „ m „ .. 

Y h(uj) + -m 2 log(Y \ a i | 2 ) + 2 m2 _ 2 m 

i =0 i =0 

if furthermore f(X) is irreducible, then Ci(u) can be replaced by 

( m—l r, m Q 1 

E h(X) + g ml °s(Y + 2 m2 “ Y 

i =0 i=0 

Finally, if f(X) is an irreducible polynomial over Q having only real 
roots, then we can choose 

Ni(u) = (2m + 1 ) 3m F(r) 4m - 1 C' 2 (w). 




Proof. Let ai,a 2 ,..., a m be the roots of f* such that |aq| > \aj\ for 
any 2 < j < m. Note that they are all distinct and also the roots of /. 
We also note that ctq is a real number. 

Denote 


r(u) = 2 _m ( m_1 h m_2 )/ 2 ( m _j_ 

Then, by (2.10) for any 2 < j < m we have 


(3.1) |aq| — | (Xj | > r(w). 

As mentioned before, for any integer n > 0, u n can be expressed as 

m 

Un = Y b i a F 
3 =1 


where b\,... ,b m are all non-zero complex numbers. Under the assump¬ 
tions, we know 


lim — = b\. 

n—>oo oOf 

Since the sequence {u n /af} is a real sequence, b\ must be a real num¬ 
ber. 

Now, we want to find a lower bound for n such that 


(3.2) l&ictfl > Y\ b i a lV 

3 =2 

So, the sign of u n is the same as that of biaf when the index n is 
greater than this lower bound. Then, everything is done. Note that it 
is equivalent to require that 

m 

N > (|oq|/M) n , 

3 =2 
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which, by (3.1), follows from the inequality 


(3.3) |6i| > ^1^1(1-r(M)/|ai|) n . 

3 = 2 

On the other hand, for any 1 < j <m, by (2.5) we know that 

|log|M<[Q(&i) :Q]h(y. 

Since bj E Q(u 0 , ..., u m _i, «i,.. •, a m ) by (2.13) and [Q(«i,..., a m ) : 
Q] < ml (ml is the factorial of m), we have [Q (bj) : Q] < ml ■ d. So 

| log |&j|| < ml ■ dh(bj). 

Using (2.11), we get 

|l°g|M < B(u), 

that is 


(3.4) exp (—B(u)) < \bj\ < exp (B(u)) 

for any 1 < j < m, where 

(3.5) 

( m —1 q m “11 

Y, h(Mj) + -m 2 \og(Y\ a *i I 2 ) + 2 m2 ~ 2 m _ 1 

i=0 i =0 

Thus, the inequality (3.3) is implied in the following inequality 
exp(— B(u)) > mexp(B(u)) (1 — r(u)/\a\\) n , 


which is equivalent to 


n > 


2 B(u) + logm 


-log(l -r(u)/|a 1 |)' 
By Lemma 2.8, it suffices to ensure that 

2 B(u) + logm 


n > 


- log(l - r(u)/( 1 + H(f *)))' 

Using the Taylor expansion — log(l — z) = z + z 2 / 2 + z 3 / 3 + • • • for 
\z\ < 1, it suffices to require that 

2B(u) + log m 
r(u)/(l + H(f*))' 

Thus, we get the following lower bound implying the inequality (3.2) 

(3.7) n > 2 m (™-l)(m-2)/2 (m + 1 )m(m-l)+6_ H -( f*ym(m-l)+2 C ^ u y 


So, if integer n satisfies the inequality (3.7), then the sign of u n is 
the same as that of bia™, especially u n ^ 0. This completes the proof 
of the first desired result. 
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When / is an irreducible polynomial over Q, one can prove the re¬ 
maining desired results by applying the same arguments and using (2.8) 
and (2.12). □ 

From the above proof it is easy to see that only with respect to SP, 
Theorem 3.1 is also true for non-real simple LRS of algebraic numbers. 
This also holds for Theorems 3.2 and 3.3. 


3.2. Characteristic polynomials with exactly two roots of max¬ 
imal modulus. The case when the sequence {w n } has exactly two 
characteristic roots with maximal modulus is quite different from the 
above case and is more complicated. 


Theorem 3.2. Let {u n }, f(X), f*(X),d be defined as in Theorem 3.1. 
Suppose that f has exactly two roots, say oq and a 2 , with maximal 
modulus such that a\/a 2 is not a root of unity. Then, for the sequence 
{u n }, SP is decidable by only testing the first [W 2 (m)J terms, where 

N 2 {u) = 2 C 3 (u) log Cfiu) 

and 

Cfiu) = F(m,d)H{f*) m3 - m2 - m / 2+3 - 5 C 1 (u), 

F{m, d) = 2 m ( m-1 )( m-2 )+ 42 7i-( m ! . df-\m + l)™ 3 /4-3m/4+6.5. 

if furthermore f(X) is irreducible, then Cfiu) can be replaced by C 2 (u). 
Here, Cfiu) and C 2 (u ) have been defined in Theorem 3.1. 


Proof. Under the assumptions, we write f(X) as 

/(X) = (X - ai )(x - a 2 ) • • • (X - a m ) G Q[X] 

such that a\,a 2 ,... ,a m are all distinct and |aq| = \a 2 \ > \aj\ for any 
3 < j < m. By (2.9), for any 3 < j < m we have 

| ay | — | oij | > s(u), 


where 

s(u) = 2( m -f- ]_^-m 3 /4+3m/4-6 +m 2 +m/2-2 

Note that for any integer n > 0, u n can be expressed as 


Un = '£ b 

3 =1 


■j Q£j , 


where b \,..., b m are all non-zero complex numbers. 

In the sequel, we want to find a lower bound for n such that 

m 

|&ictq + b 2 afi\ > J2\bj a ]l 

J—2 


(3.8) 
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So, whenever the index n is greater than this lower bound, we have 
u n 7 ^ 0. This can complete the proof. 

The key step is to get a lower bound for the left-hand side of (3.8) 
by using Baker’s inequality on linear form (2.14). Then, let the right- 
hand side of (3.8) be less than the lower bound, this can give the desired 
lower bound for the index n. 

For n > 0, we have 

(3.9) |Mi + b 2 o%\ = | 6 i<l ’ 

Here, for n > 0 we put 

A = (- i ) • r ■(-)"-1 

b i a i 

and 

(3.10) An = log(A„ + 1). 

Then, by definition, there exists an odd integer a such that |a| < n+ 1 
and 

An = a log(—1) + log(6 2 /6i) + 77 log(a 2 /«i), 

which gives a linear form in the logarithms of algebraic numbers. We 
need to handle the exceptional case when A n = 0 . 

In the following, we assume that 

(3.11) |A„| < 1/2. 


(-i)-T(-r-i 

bi 


If this is not true, then later on one can see that this implies much 
better results; see (3.12) and (3.19). 

Notice that for any complex number z with 0 < \z\ < r < 1, using 
the Taylor expansion, we have 


]og(l + z)\ = \z~£ + j~** 

2 

< ll + - + -—I— 
- 1 2 3 


i | |log(l -r)| 

\z\ =- z 


Using this estimate together with (3.10) and (3.11), we obtain 


(3.12) i|A„,| = i|log(A n + l)| < |A n |. 

Suppose that A„ = 0. Then A n = 0, that is + 62^2 = 0- Let 


A Q(ll0) • • • j 'U'm— I; ^1) • • • j ®-m) • 
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Then, [K : Q] < ml ■ d. If a\/a 2 is not a unit of K , then there exists 
a prime ideal p in the ring of integers of K such that ord p (ai/a 2 ) is 
non-zero. Then, we get 

(3.13) n < n|ord p (ai/o; 2 )| = |ord p (6 2 /&i)l < |ord p (&i)| + |ord p (6 2 )|- 
On the other hand, by definition, for any 1 < j < m, we know that 

\bj\ f = p -° rd p( fe d/^, 


where p is the underlying prime number of p, and e p is the ramification 
index of p over p. Using (2.5), we obtain 




1 log l^yl 

lpl - dp 

Wl), 

where d p = 

[K p : Qp], 

Notice 

that e p < 

: dp. So, for any 1 < j < m we 

have 





(3.14) 

|ordp(6j)| 

d P log p 


where B(u) 

has been 

defined 

in (3.5). 

Combining (3.13) with (3.14), 

we get 





(3.15) 


n < AB(u 

)• 


Now, we suppose that A n = 0 and a\/a 2 is a unit of K. Since q>i/q> 2 
is not a root of unity, there exists an embedding o : K y C such that 
|cr(ai)/cr(a 2 )| > 1. By (2.9), we have 

|cr(ai)| - |a(a 2 )| > s(u). 

Notice that 

\a(b 2 )/a(bi)\ < exp(2 B(u)), 
which can be deduced similarly as (3.4). In view of 
a(bi)a(ai) n + a(b 2 )a(a 2 ) n = 0 , 

we deduce that 

(3.16) \<j{ol\)/ cr(a 2 )\ n = |cr(6 2 )/cr(6i)| < exp(2 B(u)). 

Since 

I a( ai )/a{a 2 )\ n > (1 + s{u)/\a{a 2 )\) n > (1 + a(ti)/(l + ff(/*)))", 

where the last inequality follows from Lemma 2.8, we consider the 
inequality 

(1 + s(m)/(1 + H(f*))) n > exp(2 B(u)), 

2 B{u) 

log(l + s(n)/(l + 


which gives 


n > 
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Since log(l + z) > z — z 2 /2 > z/2 for 0 < z < 1, it suffices to require 
that 

(3.17) „ > + 

s(«) 

Notice that the lower bound in (3.17) is much larger than the upper 
bound in (3.15). Thus, if integer n satisfies (3.17), then the inequality 
in (3.16) is not true, and we must have A n ^ 0. 

Now, we assume that n satisfies (3.17). So, A n ^ 0. Applying 
Baker’s inequality (2.14), we find that 

(3.18) |A n | > exp ( —2 38 D 2 AiA 2 A 3 log (eD) log(en + e)) , 

where D is the degree of the number field generated by 62/^1 and 012/011 
over Q, and 

Ai = 7 r, 

A 2 > max{Dh(b 2 /b 1 ), | log( 6 2 /&i)|}, 

A 3 > max{Dh(a 2 /ai), | log(a 2 /ai)|}. 

Since both b 2 fb\ and ol 2 /ol\ are contained in K, we have 

D < [K : Q] < ml ■ d. 

By (2.11), we get 

Dhih/h) < D( h( 6 i) + h( 6 2 )) < 2 B(u). 

In addition, we note that 

| log( 6 2 / 6 i)| < log \b 2 /bi\ + 7 T < Dh(b 2 /bi) + n < 2 B{u) + 7 r. 

Thus, we choose 

^2 = \c,(n), 

where Ci(u) has been defined in Theorem 3.1. It is easy to see that 
A 2 > 2 B{u) + 7 r. 

Now, we want to choose A 3 . Since 

h(a 2 /a:i) < h(a 2 ) + li(«i) < 2 log M(f*) < 2 log H(f*) + log [m + 1 ), 

where the last inequality follows from (2.1). On the other hand, we 
have 

| log(a 2 /ai)| < log \ot 2 / ol\ | + 7 r < 2 log (H(f*) + 1) + n, 
where the last inequality is derived from Lemma 2.8. So, we can choose 
A 3 = m\ ■ 2 d(m + l) 0 ' 5 H(f*) 0 ' 5 . 

Define 

c = 2 39 D 2 A 1 A 2 A 3 log(eD). 
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Then, under (3.17), the inequality (3.18) becomes 

(3.19) |A n | > exp(—clogn), 
which, together with (3.9) and (3.12), implies that 

(3.20) |&ictq > | exp(—clogn). 

Substituting relevant results into the definition of c, we redefine 
c = 2 40 7r(m! • d) 3 ' 5 {m + l) a5 i/(/*) a5 C y i(u). 

Now, we are ready to find a lower bound for n such that 

m 

| Mi + M2I > ^ I Mil- 

3 =3 

This is implied in the following inequality by using (3.20) 

^ m 

-|Mrl exp(—clogn) > ^ \b] a ]\- 

3 = 3 

That is, we need that 

m 

|6i| exp(-clogn) >2^ bj\ (||/1«i|)^ , 

3 = 3 

which follows from the inequality 

m 

(3.21) |6i| exp (—clogn) > 2 E N (1 — s(n)/|ci'i|) n . 

3= 3 

By (3.4), the inequality (3.21) is implied in the following inequality 
exp(— B(u) — clogn) > 2mexp(i?(n)) (1 — s(u)/\ai\) n , 
which is equivalent to 

—nlog(l — s(n)/|ai|) — clogn > 2 B(u) + log(2m). 

As in (3.6), it suffices to reqirire that 

(3.22) ns(u)/(l + H(f*)) — clogn > 2 B(u) + log(2m). 

Notice that the integers n satisfying the following inequalities also 
satisfies (3.22), 

( clogn < ns(u)/(2 + 2H(f*)), 

(3.23) { 

I ns{u )/(2 + 2H (/*)) > 2B(u) + log(2m). 
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Since the function x/logx is strictly increasing when x > 3, for 
A > 3, if x > 2 HlogH, then x/ log a; > A. Thus, if 

(3.24) n > 4cs(it) _1 (l + //(/*)) log (2cs{u)~ 1 {l + if (/*))) , 

then the first inequality in (3.23) holds, and in fact the second one also 
holds. Note that the lower bound in (3.24) is much bigger than that in 
(3.17). 

Thus, by (3.24) we get the following lower bound for the index n 
implying the inequality (3.21) 

(3.25) n > 263 ( 11 ) log 63 ( 11 ). 

So, if integer n satisfies the inequality (3.25), then we have 

m 

|&ia"+ M 2 I >J2\ b j a ?\- 

3 = 3 

Thus, u n 7 ^ 0. This completes the proof of the first desired result. 

The remaining desired result follows from (2.12). □ 

Following the ideas in [18] and in the proof of Theorem 3.2, one can 
get an explicit result for Lemma 2.5 in the case when the sequence {u n } 
is simple and has exactly three roots with maximal modulus. However, 
this might be very complicated and deserves studying as a separate 
project. 

3.3. Quadratic characteristic polynomials. If the sequence {it n } 
is of order 2 , we can get better results. 

Theorem 3.3. Let {u n } be a real simple LRS of algebraic numbers 
defined by (1.1) of order 2. Let f (X), f* (X), d be defined as in Theorem 
3.1. If f has a dominant root, then for the sequence {u n }, SP, PP and 
UPP are all decidable by only testing the first ]_A ^3 ( 11 )J terms, where 

N 3 (u) = 4 dH(f*) 2 fh(«o) + h(in) + | log(J^ \a *| 2 ) + . 

Otherwise if f(X) does not have a dominant root and the quotient of 
its two roots is not a root of unity, then SP is decidable by only testing 
the first |_iV 4 (i/)J terms, where 

Nfiu) = 8 d ^h(ito) + 11 ( 11 !) + ^ logE l a i I 2 ) + ^ • 


21 


Proof. Let a i and a 2 be the two characteristic roots of {u n }. First, we 
suppose that |ai| > |« 2 1■ Note that in this case both ol\ and a 2 are 
real. By Lemma 2.9, we know 

|«i| - |a 2 | > 1/tf (/*)• 

Then, using Lemma 2.13 and applying the same arguments as in the 
proof of Theorem 3.1, we get that if integer n satisfies 

n > N 3 (u), 

then the sign of u n is the same as that of 61 af, especially u n 7 ^ 0. This 
completes the proof of the first part of the theorem. 

Now, we suppose that |«i| = |a 2 |. If ai/a 2 is a unit of K, then in 
view of \q>i/cx 2 \ = |a 2 /ai| = 1 , we find that a\/a 2 is a root of unity, 
which contradicts the assumption that a\/a 2 is not a root of unity. So, 
j a 2 must be not a unit of K. Note that for any integer n > 0 we 
have 

u n = biaf + b 2 a .2 

for some non-zero complex numbers b\ and b 2 , similar as the deductions 
of (3.15) we conclude that if integer n satisfies 

n > iV 4 (n), 

then u n 7 ^ 0. This completes the proof of the theorem. 

□ 

3.4. Decidability for almost all sequences. As mentioned before, 
although Theorems 3.1, 3.2 and 3.3 only concern real simple LRS of 
algebraic numbers, the arguments are still valid for arbitrary simple 
LRS of algebraic numbers when only considering SP. For example, in 
this case it is not needed that the coefficient b\ in the proof of Theorem 
3.1 is a real number. We merge them into the following. 

Theorem 3.4. Let {n n } be a simple LRS of algebraic numbers defined 
by (1.1) of order m > 2, and let f(X) be its characteristic polyno¬ 
mial. Suppose that f has a dominant root or has exactly two roots with 
maximal modulus (their quotient is not a root of unity). Then, SP is 
decidable by only testing finitely many initial terms, whose amount is 
bounded explicitly by the initial data as in Theorems 3.1, 3.2 or 3.3. 

Now, we are ready to state and prove the key result. 

Theorem 3.5. For almost all LRS of algebraic numbers with integer 
coefficients (or rational coefficients), SP is decidable by only testing 
finitely many terms, whose amount is bounded explicitly by the initial 
data as in Theorems 3.1, 3.2 or 3.3. The same is also true for both PP 
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and UPP for real LRS of algebraic numbers with integer coefficients 
(or rational coefficients). 

Proof. Let {u n } be a LRS of algebraic numbers defined by (1.1) of 
order m > 2, and let f(X) be its characteristic polynomial. 

Suppose that / is an integer polynomial. Note that / is monic. Since 
almost all monic integer polynomials are irreducible (for example, see 
[6]), we assume that / is irreducible. Then, {u n } is a simple LRS. 
By Lemma 2.3, almost all monic integer polynomials have a dominant 
root. Then, this case of SP is done by applying Theorem 3.4. 

Now, assume that / is a polynomial over Q. Let S be the smallest 
positive integer such that 5f(X ) G 7L[X}. 5f has the same roots as /, 
but note that 5f may be not monic. Since almost all integer polyno¬ 
mials are irreducible (for example, see [15]), by Lemma 2.4, we only 
need to consider the case when Sf is irreducible and non-degenerate. 
Furthermore, by Lemma 2.4, we can assume that 6 f has a dominant 
root or has exactly two roots with maximal modulus. Then, this case 
of SP is also established by using Theorem 3.4. Thus, this completes 
the proof of the first part. 

For real LRS of algebraic numbers with integer coefficients (or ratio¬ 
nal coefficients), the remaining desired results can be proved similarly 
by using Lemma 2.6 and Theorems 3.1 and 3.3. Note that for an 
irreducible and non-degenerate characteristic polynomial /, if it has 
exactly two roots with maximal modulus, then these two root must be 
non-real, thus this case has been done in Lemma 2.6. □ 

3.5. Estimate the exceptional set. Regarding the results in Theo¬ 
rem 3.4, one question of interest is for how many such LRS the Skolem 
Problem might be not decidable explicitly. It depends on the way how 
to count LRS. 

For LRS of algebraic numbers with integer coefficients of order m > 
3, we define that two LRS are equivalent if they have the same char¬ 
acteristic polynomial. So, it is equivalent to ask for how many such 
equivalence classes the Skolem Problem might be not decidable ex¬ 
plicitly. For the equivalence classes corresponding to the set S rn {H) 
(defined in Section 2.1), let A m (k,H) be the number of polynomials 
/ G S m (H ) having exactly k roots with maximal modulus, then by [9, 
Theorem 3.2] and [7, Theorem 1.1] we have 

(3.26) \S m (H)\ - A m ( 1 , H) - A m {2, H) < Cl (m)H m -\ 

where C\(m) is a function with respect to m. According to Theorem 
3.4, this can be viewed as an answer to the question. 
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On the other hand, for LRS of algebraic numbers with rational co¬ 
efficients of order m > 3, we can define the same equivalent relation. 
For the equivalence classes corresponding to the set 

let A^(k, H ) be the number of polynomials / G T^(H) having exactly k 
roots with maximal modulus, then by [9, Theorem 3.4] and [7, Theorem 
1.4] we have 

(3.27) \rjH) | - - Al(2,H) < 

where C 2 (jn) is a function with respect to m. This can be viewed as 
another answer to the question. 
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